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ABSTRACT 



For many years the cylindrical dipole and circular loop have been 
thought of as the basic radiating elements. Due to its relative 
simplicity the study of the cylindrical dipole is quite complete, but 
owning to the degree of complexity, progress in the theory of the loop 
has been much slower. In recent years the isolated circular loop 
antenna of arbitrary size has been the subject of mathematical and 
experimental analysis [7], [6]. It was not until 1970 that the first 
complete and general mathematical solution for the mutual admittance 
of two arbitrary size coaxial circular loop antennas was published. 

The problem was formulated by Adler in terms of coupled integral equa- 
tions for the current distributions on the loops [1]. The equations 
were solved using the Fourier series representation of the current 
distributions. The study presented here is the next logical step 
beyond the work of Adler. Using his current equations the current 
distributions are calculated at each point around the two loops. The 
input admittance, electric field and power gains are calculated. 
Radiation patterns and the loop current distributions are plotted for 
two cases: (1) both loops driven and (2) one loop acting as a short 
circuited parasite. 
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I. INTRODUCTION 



A. HISTORY OF WORK ON THE LOOP ANTENNA PROBLEM. 

In 1897 Pocklington excited a thin loop by a plane wave and 
obtained a solution for the induced current in the form of a Fourier 
series expansion [9]. Hallen continued the work by producing an 
expression for the current distribution in the driven case [3]. Due 
to singularities arising in the Fourier series coefficients, Hallen's 
technique applied only to electrically small loops. Storer attempted 
to surmount the limitations of Hallen's work and derive an expression 
for the current in the case of a large loop [10]. He expressed the 
series in integral form and evaluated the integral in the sense of the 
Cauchy principal value. Experiments of Kennedy showed Storer' s cal- 
culations produced an unacceptable deviation between calculated and 
measured susceptance [6]. Wu reviewed the work of Hallen and Storer 
and discovered that the shortcomings arose from over-simplifying the 
problem [11]. That is they assumed the electric field intensity on 
the surface of the conductor was uniform around the circumference with a 
filamentary current flowing at the center of the conductor. 

Wu then modified Storer's solution by changing the one dimensional 
integral equations to two dimensional integral equations in order to 
account for the electric field variation around the circumference of 
the conductor. Wu assumed a longitudinal current flow on the surface 
of the conductor and allowed no circumferential component. He showed 
this to be equivalent to a uniformly distributed current density on the 
cross section of the conductor, which allowed the current to be treated 
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as filamentary and concentrated at the center of the conductor. After 
this simplification the equations were again one dimensional but no 
longer required uniform electric field around the circumference. King 
applied Wu's modification to calculate the admittance of a single loop 
[7]. This proved Storer's work to be valid but only for the case of 
small loops. 

The first attempt at solving the multi-loop problem was by Iizuka 
for the special case of two identical loops [5]. In 1970 Adler advanced 
King's work by devising a completely general mathematical solution to 
the mutual admittance of large arbitrary size coaxial loops. Adler also 
measured the admittance in the two-loop case and found close agreement 
with theoretical values [1]. 

B. SCOPE AND LIMITATIONS OF THE STUDY 

This paper completes the study of large arbitrary size coaxial loop 
antenna arrays. Current distributions and radiation patterns are 
calculated and plotted for various array geometries. Although only the 
two-loop situation is presented herein, the theory applies to the 
n-loop situation as well. Since the $ (polarized) mode electric field 
component is of primary interest, this mode is examined in detail. The 
e mode is calculated and one representative field pattern is plotted. 

The range of loop dimensions selected is for the resonant and near 
resonant case, a, b, and c are the conductor radius, loop radius and 
loop spacing respectively; the range of parameters are expressed as 

0.5 s kb i 2.0, 0.5 i kc i 7.0 

where: 
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ka is determined by the thickness parameter: 

n = 2 In 2llb and \ is the wave length, 
a 

Having 7.5 s n s; 15 represents a ratio of b/a from approximately 10 to 
300. The above parameter limitations are a result of the approximations 
used in the theoretical analysis of the Fourier coefficients of the 
current distributions. 

Although perfect conductors in free space are specified, the 
effect of lossy media can be included merely by considering the propaga- 
tion constant complex. The loops are driven by delta function voltage 
generators, symbolized by V 6(<t>). The two coupled antennas are treated 
as a two-port network and conventional circuit analysis techniques are 
employed to produce a convenient matrix of equations. 
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II. THEORY 



The conductor radius is assumed to be small relative to the loop 
radius, which allows the current density in each loop, I-j (<t> ) and 
to be approximated by an equivalent one-dimensional filamentary current 
distribution at the center of the loop. The electric field is derived 
by the scalar and vector potential in terms of the distribution of the 
current. Readers are referred to Adler's work for the mathematical 
details of the current equation formulation. Only an outline of the 
solution is presented here. 

As shown in Figure 1, the cylindrical coordinates for points on 
loops 1 and 2 are (b-j , 4> ' ,o ) and (b 2 ,<J>,c) respectively. The angle <f>' 
represents the coordinate of the current element, while the angle <f> 
represents the electric field observation point coordinate. The 
distance from the source current element on one loop to the electric 
field observation point on the other loop is: 

Rl2 = R21 = V(bl" b 2) 2 + + 4b-jb2Sin 2 • 

The distance from one point to another on loop 1 is: 

R-|-| “V" 4b^ sin 2 + 4a^ sin 2 

as shown in Figure 2. The same expression applies to loop 2 if the 
subscript 1 is changed to 2. ip is the angular distance around the cir- 
cumference of the conductor surface. In order to simply the notation 
in this development, the notation for harmonic time variation, e Ja)t 
has been suppressed. 
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Figure 1. Geometry of Two Unequal Loops 
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Figure 2. Geometry of Points on a Loop 
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e is the total electric field at a point due to loop currents 
I-j ( 4> ) and I 2 (<|>). The boundary conditions on e at the conductor 
surface require the tangential component to be zero at the surface; 
since no circumferential current is assumed, the E^ component is the 
only tangential component present. The delta function driving 
generator placed in the loop causes a singularity in the electric field 
at the gap such that: 



and hence: 




(1) 



The boundary conditions on both loops are analogous, is the <j> 
component of the total electric field due to the currents on both loops. 
These boundary conditions, when applied to the expressions for the total 
electric field intensity, specify the following coupled complex integral 
equations for the unknown current distributions: 



V ] 6(*) = + K 12 (4>-* , )I 2 U , )]d*' 



*n Z * 11 

v 2 6(«J>) = ^2. J [K 21 (*-* , )I 1 (+') + K^U-m^'nd*' 



-n 



In these equations 



KijU-*') = [kb.jbj cos(<f>-<|> ' ) + £■ 



( 2 ) 

( 3 ) 

( 4 ) 
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where K. . is the kernel of the integral equations and 

' J 

[kb,b. cos(<j>-<j>') + i iL. ] 

J K 3<t» 2 

operates on W.j (<{>-<}>'). Wu's contribution to the solution resulted in 
integrating the term: 



j-jkR-j j 



R ii 



to 



n -JkR f1 



/■« -JKK 

get: «„<*-*') =-_y — 



dip 



(5) 



where W^ (<}>-<{>') is the Fourier coefficient of: 



R n 



When i f j 



-jkR. . 

W . . (<i>-4> ' ) = - U- • 

1J R. . 

1 J 



( 6 ) 



n 0 is the characteristic impedence of free space. 

Equations (2) and (3) are now solved for I-jU) and ( 4>) by the 
orthogonal expansion technique. The kernels of the integral equations 
and the currents are specified in terms of their respective Fourier 
series as follows: 

°° . -jn<j>' 

1.4) - l r e (7) 

n=-co 



K. ) = l 

n=-°° 



.. - jn 1 ) 

V e 



( 8 ) 



By applying equations (7) and (8) to equations (2) and (3), inter- 
changing the order of integration and summation and performing the 
integration over d>, equations (9) and (10) 'are formed: . . 
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( 9 ) 



11 1 1? o j n 'f> 

* t” 1 to 1 )e •»,«(♦) 

2 n n n n n 1 



j2» J [o 21 ,1 + a 22 ,2 ] e J "* . v s(+) 

2 p; n n n n 2 y 



( 10 ) 



After multiplying both equations by e” jn ^ and integrating over <j> from 
-II to n, the following simultaneous equations are obtained: 



11 


12 


n 


a n 


21 


22 


n 


cu 

n 



_ 1 



jnno 



'1 



( 11 ) 



1 o 

Equations (11) are solved for 1^ and I n and substituted back into 
equation (7) to get: 



and 



1.4) = -J— Icos n<ji 
1 jn no n 



i 2 (*> 






£cos n<j> 





12 

a n 


V 2 


a 22 

n 


11 


12 


(X 

n 


a n 


21 

a n 




•i 1 


V 1 




V 2 


11 

a n 


J 2 

n 


a? 1 

n 


a 22 

a n 



-jn* 






( 12 ) 



(13) 



The task now is to evaluate the a^'s. The coefficients a]3 are put in 

n n 

terms of the Fourier coefficients k 1 '3 of W* (<f>-<|> 1 ) by letting: 

n J 
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u (A a'\ - v 1 J ' ) 

U-<f> ) = I < J e 

J n=-» n 



( 14 ) 



where: 



^ = k’J = r W, .(*-*') eJ n (+-*') d*. 
n -n 2lt J ^ 



(15) 



-n 



The following difference relationship is arrived at after differentia- 
tion of equation (15): 



ij _ ij _ kb 1 *b i 

a - a - — 1 — J- 

n -n 2 



k 1 ?. + ^ J i 

n+1 n-1 



JL_ 

k n 



. (16) 



a = kb.b. K-, 
o i J 1 



1J 



(17) 



For the evaluation of refer to appendix A. 



The coordinate system is as indicated in Figure 3. The expressions 
for the far-zone electric field and gain are developed around Adler's 
coupled current equations (12) and (13) rewritten below: 



jn no 



V a 12 
1 a o 

22 



v i “J 2 

22 



\1 a dd “ V a“ 

11 22 12 21 + JjL cos(n<|,) II 22' a l2 a2l 

a a -cu an- a a*- 4 - -a n a{r 1 

o o ° o n n n n j 



(18) 



i 2 (*) = -4 



jnn 0 



■l 1 » 

a? 1 V 



aj 1 a 22 -aj 2 a 21 n=l 



+ cos(n<j>) 



ajl V , 

a 21 Vo 
n ^ 

a 11 a 22 -a 12 a 2 1 
a n a n a n a n j 



.(19) 



Let I-| (<}>), the current in loop 1 conductor element b df at the angle <j>' 
be as given in (18) above and ^U) be the same except for loop 2. 
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Figure 3. Vector potential for one loop due to element of 

current I(<j>') b d<j>' (a) at arbitrary point and (b) 
at point in plane of loop; (c) spherical coordinates. 
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The element of vector potential dA^ at the point (x,y, 
1 1 ( 4> * ) is composed of two components, dA-j x and dA^ where: 



and 



dA-j x = -dA-j sin <j> ' 



dAjy = +dAj cos <j)' 



So that: 



-jkR 



dA-j x = - 1L_ I-, (4> ' ) sin <f,' b-jd^, 1 



and 



-jkR 



dA ly = + 4n M*'* ^“r — cos b l d +' 



where 



= Vz^ + s 2 = V (r cos e) 2 + s 2 



and 



s^ = (r sin e)^ + - 2br Q sin e cos (<j>-<j>') 



o 2 

When b «r Q » R - r 0 -b sin e cos ' ) . 

2 2 

The vector potential from loop 1 at points where b <<r 0 is 



by: 



a r _ _ b l y e"^ kr o f 2 j / ,x jkb, sin 0 cost*-*') 
A lx ~4n r 0 J ¥♦ } e 1 

/ 2n 

Ij(*') e jkb l sin 6 cos( *“* ,) c 



sin 



A r.^e^ r ' 

ly 4n r „ 0 



The same procedure is followed to get the vector potential 
by loop 2. 



z) due to 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

thus given 

<j> , d<j>' (26) 

$'d$' . (27) 
radiated 
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Equation (18) is substituted into equations (26) and (27). The 



resulting integrals are solved as follows: 



Ii (o) = - 



nc 



v i »J 2 

V 2 “n 2 



*— n — * 



(28) 



where I-j (o) is the nth term in the series for loop 1 current at <j>=o 
n . 

(the feed point) and c = 4i/e where e is permittivity. 



and 



11 22 12 21 
“n “n a n a n 


(29) 


b l u P -Jkr_ 

= Tf -T- 3 - [ 1 (1) ' 
HJ1 r o n 


(30) 


kb-j sin e . 


(31) 


and costj) 1 in exponential form so that 




Qi 

= 3f (-k. * V 

1 


(32) 


Ql 

= 4n (K + K b ) 

j2 a D 


(33) 



where 



2n j(n+l)<j>' jB cos ( ' ) 

K a = / e e 1 d<j>' 



Kl 



■/ 

o 

JZu j(n-l )<j> 1 jB-j cos ( 4>-4» ' ) ^ 



d<j) 1 . 



(34) 



(35) 



These integrals are of the form of the Sommerfeld representation of the 
Bessel function: 



j-n m j 

V B > - hrj e 



jBcos$> jn$ 

e d$> 



(36) 
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Letting $ = <p-tp' gives: 



.-n /-2n-<j) jBcos ( <^~ 4 > ’ ) jn^' 

° n e e d *' 






(37) 



This requires that 



K g = 2 ji j 



n+1 j(n+l)<j> (B-i ) 

e J 



n +1 



and 



K h = 2n j 



n -1 j(n-l )<j) 



(B-i ) 



J n -1 



(♦) 



(38) 



(39) 



The components of the vector potential due to the component I-j in 
spherical coordinates are: 



and 



AC = (Af cos d> + At" sin d>) cos e 
\ x n 'y n 



AT = (-AC sin <j> + A] cos <|> ) . 
4 > n x n Yn 



Substituting (30), (34), and (35) into (40) and (41) gives: 

Ml. 



A 



e n 2 



— [(Kg -K b ) cos<() - j ( Kg +K b )sin<j)] cose 



J’Ql 



11 [K a e" J< ^ -^e^] cos 4 ) 



Af = 5ln [(Kg +K b ) cos<j, - j ( K a -K b ) sin<t.] 
9n 2 



Qi 



2 (K a e" J<!) +K b e J ^) . 



Substitute (38) and (39) into (43) and (45) to get: 

A le ’ Q 'n ej " U+ " / 2 )[J n+l (Bl) + J n- 1 <B,) ] cose 



= -2n Q, e j n( * +II/2) 2_ J (B 1> 



cose 



(40) 



(41) 



(42) 

(43) 

(44) 

(45) 



(46) 

(47) 



and 
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AV (B,). ( B)] 

<t>n n 

A^ = -j2n Q 1n e J ’ n(<j>+n/2) J^ B 1^ . 

By summing all of the contributing vector potential components and 
using the recurrence relationships: 



(48) 

(49) 



2jl j ( x ) = j ( x ) + j (x) 
x n n+1 n-1 



and 



2J ' 



,(x) =J (x) 



n-1 



- J 



(x) 



n+1 



the vector potential at a point due to loop 1 is: 

-k =■ 5 \.‘""+fc V*,‘" 

0 -co 0 -co n C 1 



cose 



(50) 



and 



r 00 „ 00 jnU+n/2) (b \ 

A 1 -I*U - -JZH Z <h n e J^l' . 



<t> -» y n 



b -jkr 0 

Substitute Q-| = -=3ii tH into equations (50) and (51) 



4n 2 c 



A" - -2n Q, l 

0 -oo 



a!| = -j2n Q-, l 

(J) -oo 



v i -l 2 

h % 



jn( 4 >+IL) J ( kb, sine) 
, 2 n n i 



kb 



cose 



1 sine 



n 

h 4 2 



(kb-j sine) 



(51) 



(52) 



(53) 



Since = a^ 1 and J_ n ^ x ^ = (-l) n J n ^ X ^> equations (52) and (53) can 
be rewritten as: 
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C - 



-2n Q, J n 

1 



V 1 “J 2 

v 2 4 * 



O n (kb 1 sine) 



cose 



D kb i sine 



jn(<f>+n/2) (-l) n - j n ( <j>+n/ 2 ) 
[e e ] 



( 54 ) 



and 



A r 

A 1 



* 



= -J‘2n Q, 



vi 2 

V 2 4 2 



^ 0 ( kb-j sine) + l 

1 



v i 

v 2 a n 2 
D_ 



( kb-j sine) 



jnU+n/2) -jn(<j>+n/2 

[e +(-!)" e ] 



(55) 



The electromagnetic field from loop 1 is then given by: 

e!f = -jto(e Af + l A[ ) . (56) 

King shows that for a single loop of moderate size, | kb | <2 . 5 , the 
significant contributions are from the first three terms of the series 
[7]. 

The electric field patterns for the <j> and e modes for loop 1 are 
readily obtained using the first three terms of the respective series. 
The identities: 

and 

J,(X) . J (X) . J, (*> 

X 

are used. 
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n=0: 



A l e = 0; A l, = J2n Q-j 
e o *o 



Vi *J 2 

v 2 a o 2 
D_ 



J.| (kb^ine) 



£ lo = * \ s2i v 



V 1 “o 2 

£ 



V2 <*22 



Jl ( kb-j s i ne ) 



D, 



(Note that for the constant current distribution situation the 
has only the n=0 term.) 
n=l : 



A v = 4n Q 1 



1 



v 12 
V 1 a l 



/ 22 
'2 a l 



J 1 (kb-|Sin9) cos0 S1 - n< j ) 



kb-! sine 



%, = 4n Q i 



V, 

vf 

D, 



(kb^sine) cos<j> 



e i „ = 0 Ei +$E 



1* 



1 



E L = ' ;i411 Q i 



(jl) 



v, 4 * 

V 2 af 



J l( kb l s1 ’ nQ ) cose sin<|> 



kb-| sine 



Ei = -j4n Qia) 

91 



v i -1 



12 



22 



V 2 a l 



D 



1 



n=2: 



A l, 2 = j 4 " «1 



v,«i 2 

V 2 °2 2 
Do 



(kb-| sine) cos<j> 



2j 2 <kb lS , ne ) cose sin2+ 

kb-| sine 



(57) 

(58) 
series 

(59) 

(60) 
(61) 
(62) 

(63) 

(64) 
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(65) 



A r 



j2<f> 



= j4n Q, e 

♦2 



V, ^ 

h 4 2 

Do 



J 2 ( kb 1 sine) cos2<{) 



el - 0 Ei + <j) Ei . 

1 2 1 e 2 l(f, 2 



( 66 ) 



Ef = (4n Q^) 



V, „12 

V 2 c.^2 

D, 



2J 2 ( kb-| sine) cos6 sln2( , 
kb-j sine 



(67) 



4>2 



= (4n Q lU ) 



V 1 «2 2 

V 2 a 2 2 

Or, 



J 2 (kb^ sine) cos2<j> 



( 68 ) 



The field from loop 2 is found in the same manner. To find the 
total field at a point for each mode, each term of the series for loop 2 
for both modes is multiplied by the phase factor, exp (jkc cose) to 
account for the separation of the loops. The total electric fields for 
the $ mode and e mode are expressed as follows: 



e e = E f + E - e 



■r 

: 2a 



jkc cose 



+ Ei . + Eo e 



jkc cose 



and 



_ 1 . T to t 

i > 0 H 0 <h N>i 

r jkc cose 

+ e[ + Eo e 

<f>2 ^2 



„ „ jkc cose r „ jkc cose 

r r* * r ' , r* ' ■ r* ' 

£ e = E l 0 E 2 e e + E le, + E 2 01 e 

w o e o 1 9 1 



+E ln + E 2 e 
'02 ^©2 



jkc cose 
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The gain of the array is calculated by dividing the square of the field 
by the power radiated. Where: 

Power Radiated = P rad = v| • Y IN + VJj • Yj^ . (69) 



In order to find the input admittance, it is necessary to first derive 
expressions for self and mutual admittance. By treating the two coupled 
loops like the two-port netv/ork in Figure 4, equations (70) and (71) can 
be written directly: 



1,(0) . Y„ v, + Y, 2 V 2 


(70) 


i 2 (o) = y 2 , v, ♦ y 22 v 2 . 


(71) 




+ 



Figure 4. Two-Port Network 



I 1 (o ) and ^(o) are the feedpoint currents. Because the network is 
reciprocal, = Y 21 - Solving for mutual admittance: 



^( 0 ) 



'12 



1 



V-|=o JHn 0 



-° J2 + 2 1 -J? 



0 r 



n=l ^n 



Self-admittance is by definition: 

r il 



I,(o) 

Y, = 1 



V 1 

I 2 (o) 



'22 



v 2 =o 



V-| =0 



v 2=l 



(72) 



(73) 



(74) 
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The input admittance is by definition: 




( 75 ) 




( 76 ) 



In order to simplify the digital computer program used later in solving 
for the electric field and gain, it will be assumed that both driving 
voltages are purely real. This will allow to be replaced by Gj^ 
where equals the real part of Yj^. 

Then: 




( 77 ) 



and 




( 78 ) 



Final ly. 




2 



( 79 ) 



and: 



Gain 



e 




( 80 ) 



r 1 r r r 

Ej> and Eo represent Eo and E 0 respectively multiplied by the phase 




+jkc cose 



factor (e ) which accounts for effects of loop separation. 
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III. DISCUSSION OF NUMERICAL TECHNIQUES AND COMPUTER PROGRAM 



The terms were computed on an IBM System/360 Model 67 digital 
computer in FORTRAN IV language using Adler's basic program [1]. Beyond 
this point the programming was straight forward. Flow charts on 
Figures 5a, 5b, and 5c show the major steps in the program. The flow 
chart of Figures 5a and 5b is for computing I-|(<j>), ^U), Vi’ y s2» y ii> 
Y 22 > and where Y s -| and Y $ £ are the self admittance of loops 1 and 2 
respectively with mutual coupling effects ignored. Figures 5a and 5c 
show the steps in computing Gain^ and Gain 0 . The computer results for 
current in the isolated loop were compared with King's [8], while the 
gain was compared with the work of Harrington and Mautz [4]. In both 
cases the agreement was exact. The incrementing of parameters b and c 
in the program allows for changing loop size and loop spacing. 
Incrementing 6 and <j> allow for computation of gain patterns. 
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Figure 5a. Flow Chart of Computer Program 
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Figure 5b. Flow Chart of Computer Program 



28 



* 















Figure 5b. Flow Chart of Computer Program (Continued) 
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Figure 5c. Flow Chart of Computer Program 
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Figure 5c. Flow Chart of Computer Program (Continued) 
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IV. DISCUSSION OF RESULTS 



A. CURRENT DISTRIBUTION AND GAIN PATTERN CORRELATION 

Figures 6 and 7 show the patterns of Gain^ in the x-z plane and the 
x-y plane respectively. In all cases presented kb-|=l and loop 1 is in 
the x-y plane. Loop 2 is positioned coaxially above loop 1 at a 
distance c. Figures 8, 9, 10, 11 and 12 show the corresponding current 
distributions for the case of both loops driven in phase by one volt 
delta function generators. 

For .5 < kc < 1 and .5 < kb£ <1.5 the radiation patterns do not 
change significantly and only a slight change in phase of loop 2 current 
is noted. When kc > 1 radical changes in pattern shapes occur. This 
shows that as the separation becomes larger than about X/6, the effect 
of the array factor predominates. When kb£ = .5 the circumference of 
loop 2 is x/2. In this case the loop supports a 1/2 cycle standing 
current wave. 

Figures 13 and 14 are the Gain^ patterns in the x-z and x-y planes 
respectively, for the case when loop 2 is parasitic. Figures 15, 16, 
17, and 18 are the corresponding plots of the current distributions. 

When the two loops are closely coupled, kc = .5 and resonant the 
currents are large and almost equal in magnitude but 180 degrees out of 
phase. Even though Gain,}, is relatively high in this configuration, the 
benefit of large currents is not realized because of the phase differ- 
ence between the two currents. As the separation of the loops is 
increased from approximately x/13 to approximately X/6, the currents 
drop by a factor of about five; however, the Gain^ drops only slightly. 
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Figure 6. Gain <{> patterns in x-z plane (<j>= o°,180°); feed point 

at <}>= 0 o , e=90°; <j>=o° at right, e=o° at top; fi=10, kb-|=l, 
Vi=V 2 =l; kc=.5,l ,3,6.3; kb 2 = .5, 1,1. 5, 2. 
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Figure 7. Gain <j> patterns in X-Y plane, e=90°; feed point at 
bottom (<j>=o°); n=10, kb-. =1 , V-i=Vp=l; kc = .5,1,3, 
6.3; kb 2 = .5, 1,1. 5, 2. 
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fi=15, V,=V ? =1, kb =1 , kc=.5; (a) kb ? =.5; (b) kb 
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Figure 13. Gain <j> patterns in X-Z plane ( 4 > =o° , 1 80° ) ; feed point 
at <j>=o°, e=90°, <j>=o° at right, e=o° at top; n=10, 
kb-j =1 , V-j =1 , V 2 =0; kc=. 5,1 ,3,6.3; kb 2 =.5,l ,1 .5,2. 
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Figure 14. Gain $ patterns in X-Y plane, e=90°; feed 
point at bottom U=o°; n=10, Vt= 1, V 2 =o, 
kb-j =1 ; kc=.5,l ,3,6.3; kb 2 = . 5, 1 ,1.5,2. 
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Figure 15. Magnitude and phase for loops 1 and 2 current, distributions for 
n=15, V,=l , V 9 =0, kb =1, kc=. 5; (a) kb 2 =.5; (b) kb ? =l; (c) kb ? = 
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Figure 16. Magnitude and phase for loops 1 and 2 current distributions for 
fi=15, V,-!, V =0, kb =1, kc=l ; (a) kb?=.5; (b) kb 9 =l ; (c) kb 0 =l 
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This observation is explained as follows: although the currents are much 
smaller when c - x/6, the phases of the two currents are no longer 
exactly opposing. This in turn allows the total vector potential of the 
<J> mode in both cases to be almost equal. 



B. CURRENT DISTRIBUTIONS 



Referring again to the current plots, the following observations 
are made. When both loops are resonant and driven, the current peak 
increases as c increases from A/1 2 to A and decreases as c increases 
beyond A. When c is increased beyond 2A (not shown) the mutual coupling 
is negligible and the loop currents are very close to those of isolated 
loops. This phenomenon is due to the behavior of the input admittance, 
where: 



IN ] = Y 11 + Y 12 



V 1 = v 2 = 1. 



By computing Y ll’ Y 22’ Y 1 2 and Yj^ ^ the following was noted: (a) with 
close coupling has a large negative real part which approaches zero 

as c is increased; (b) the real part of Y^ is always positive; (c) 

Yjhj increases as c increases. Therefore, since 1(0) = Yj^V, the current 
increases as Y j ^ increases. 

Of particular interest is the behavior of the current in the 
parasitic loop. In all cases both resonant and non-resonant the 
magnitude of the parasitic current, although reduced, paralleled the 
magnitude of the current on the driven resonant loop. In the resonant 
case when the spacing was close, the two currents were 180 degrees out 
of phase, but as the spacing approached A, the currents approached the 
in-phase condition. Also in the parasitic case the phases of both 
currents tend to be constant with an abrupt reversal at the nulls. 
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C. GAIN PATTERNS 



Figures 6 and 7, the Gain^ patterns in the x-z plane and x-y plane 
respectively are again referred to. In these figures n=10, Vi = V£= 1 volt, 
kb-|=l, kc=.5, 1, 3 and 6.3 from left to right, while kb2=.5, 1, 1.5 and 
2 from top to bottom. In all patterns presented "+" marks the antenna 
center and marks unity gain. With kb 2 fixed, the patterns are 

seen to change shape as kc increases from .5 to 6.3. This is due in 
small measure to the change in mutual coupling and primarily to the 
array factor. The smallest Gain^ in the z direction occurs at kb-j = 
kb£ = 1 and kc = 3 (c * x/2); the largest Gain^ in the z direction 
occurs also with the resonant array but with kc = 6.3 (c = x). In the 
former case the vector potentials were nearly 180 degrees out of phase 
and in the latter both vector potentials were in phase. In the extreme 
case where kb 2 = 2 and kc = 6.3, lobing occurs due to the combination of 
the array factor and the diverse current distributions on the two loops. 

Gain^ for the parasitic case is illustrated by Figures 13 and 14. 

The higher gains and more interesting patterns are attained in this 
configuration. As would be expected, less lobing now occurs and the 
parasitic element functions as a reflector. Of special interest are the 
kb 2 =1.5, kc = 1 and kb 2 = 2, kc = 6.3 patterns. In the first of 
these particular arrangements the loop array pattern resembles that of a 
log periodic antenna. Figure 21 shows the effect of changing kb£ to 1.2 
and varying kc from .6 to 1.2. When kb 2 = kc = 1.2 the antenna is uni- 
directional. In contrast to this configuration, note the resulting 
pattern when kb 2 ' = 2 and kc = 6.3. The pattern is almost square with 
the gain across the top nearly uniform. Figure 14 shows a section of 
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these same patterns in the X-Y plane. These patterns closely resemble 
those of the familiar single-loop pattern in shape. 

The effect of changing the b/a ratio is illustrated in Figure 19. 
The change in Si is from 7.5 to 15, which represents a b/a ratio change 
of about 10 to 300. Although the effect is measurable, it is concluded 
that the loop array is not sensitive to changes in SI. 

Figure 20 shows the consequence of ignoring the mutual coupling 
when calculating the gain pattern (the dotted curves are with coupling 
ignored). In Figure 20a (kc = .5 and kb-| = kb 2 = 1) the coupling is 
very significant, but as seen in Figures 20b and 20c (kc = 6.3, kb 2 = 1 
and 1.5 respectively) the coupling has little effect. 

Sample patterns of Gain 0 are shown in Figure 22. Although 
significant, this mode is certainly of secondary importance. 
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Figure 19. Gain <f> patterns in X-Z plane (<j)=o o ,180°) for 
V 1 =V 2 =kb 1 =l, kc=6.3; kbo=.5,l ,1 .5,2; n=7.5,15; 
<j>=o° at right, e=o° at top. 
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Figure 20. Gain $ patterns in X-Z plane with (solid lines) 
and without (dotted lines) mutual coupling 
effects. ft=10, kb-|=l, 4>=o° at right, e=o° at 
top. 
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Figure 21. Gain 4> patterns in X-Z plane (<#> =o° , 1 80° ) for 

=10, V] = l , V 2 =o, kb-j =1 , kb 2 =1.2, kc=.6, .8,1 .2. 




Figure 22. Gain e patterns for f2=10, kb-.=l in (a) <{>=45 0 , 
225° plane; (b) <j»=90°, 270°; (c) <f>=90°, 270°. 
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V. SUGGESTIONS FOR FURTHER STUDY 



Since this study is of a theoretical nature an immediate require- 
ment is to confirm the results by experiment. First the current distri- 
butions and next the field intensity patterns for various arrangements 
should be measured. Also of interest in the two loop array is the 
consequence of open circuiting the parasitic loop. 

The current and gain equations then should be solved for the three 
element case. Of interest would be the consequence of arranging three 
or more elements to make a log periodic loop antenna. 
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APPENDIX A 



Evaluation of Fourier Coefficients 

n 

Case I: i = j 

Combining equations (5) and (15) results in: 



K 



ii 



n 




dijj de 



(81) 



By substituting A.- = 2a^ sin|- in equation (81) and interchanging the 
order of integration one arrives at the expression: 




where: 



W 




de; n > o 



(82) 



(83) 



Storer showed that for the case of A^ <<b^ : 
, /8b, \ ? n-1 1 

b iW -HxrH J, ^ 



1 

2 



2kbi 

J Cfi 2n 



(x) + j J 2 n( x ) ^ dx 



(84) 



where fl m (x) is the mth order Lommel -Weber function and J m (x) is the mth 
order cylindrical Bessel function of the first kind. 
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Wu devised an approximation method to evaluate equation (84) based 



on the facts that (a) the integral representation of: 



1 K„ 



£- 



o 1 

k 2 )2 



where K q (x) is the modified Bessel function of the second kind of zero 
order and the integral representation of b^M n (A 1 -) are similar and (b) for 
A.j small the difference between b.jM n (A^) and the small argument approxi- 
mation for K q is independent of A^ . The consequence of this is for 
moderate size n, 



b iV A i) = 



. 1 



2n 




+ 



r 2Kb,- 

2 c n J [ a 2n( 



x) + j J 2n (x) 3 dx 



(85) 



n-1 



ii- 1 -I 

where C = y - 2 Y „ , + ln(4n) and y is Euler's constant, 
n ' _e_ 2m+l 



m=o 



Because A^ is small, the K q term in equation (85) can be slightly 
modified with a negligible loss in accuracy to: 




By employing equations (85) and (86) equation (82) can be integrated 
to get: - 
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and: 



kC 11 - 1 

o 2nb i 





2kb ) 

Cfl 0 (x) = J J 2n( x ^ dx [ ' 



Case II: i^j 

By combining equation (15) with: 



R ij = [4b i b j sin 2 + (b i 




+ c 2 ]2 



( 88 ) 



and 



-jkR,-,- 
. \ _ e iJ 



w ij (♦-*') - R 



ij 



equation (89) follows: 



Jj _ 1 

n 2n 




-n 



e jn<j> e -jk [4b i b J . sin 2 | + (b i -b j ) 2 + c 2 ]2 
[4b i b J - sin 2 i + (bi-bj) 2 + c 2 ]l 



dij) . 



(89) 



It is necessary to solve equation (89) by numerical methods and the 
reader is referred to Adler's original work for the details of this 
process [1]. 
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